arXiv:astro-ph/0112348vl 14 Dec 2001 


THE PRIMORDIAL HELIUM-4 ABUNDANCE DETERMINATION: 
SYSTEMATIC EFFECTS 

TRINH XUAN THUAN 

Dept, of Astronomy, University of Virginia, Charlottesville VA 22903 

YURI. I. IZOTOV 

Main Astronomical Observatory, Kyiv 03680, Ukraine 


Abstract. 

By extrapolating to O/H = N/H = 0 the empirical correlations Y-O/H and Y-N/H defined by 
a relatively large sample of ~ 45 Blue Compact Dwarfs (BCDs), we have obtained a primordial 
4 Helium mass fraction Y p = 0.2443±0.0015 with dY/dZ = 2.4±1.0. This result is in excellent agree¬ 
ment with the average Y p = 0.2452T0.0015 determined in the two most metal-deficient BCDs known, 
I Zw 18 (Zq/ 50) and SBS 0335-052 (Zq/ 41), where the correction for He production is smallest. 
The quoted error (Ict) of <J 1% is statistical and does not include systematic effects. We examine 
various systematic effects including collisional excitation of Hydrogen lines, ionization structure and 
temperature fluctuation effects, and underlying stellar He I absorption, and conclude that combin¬ 
ing all systematic effects, our Y p may be underestimated by ~ 2^4%. Taken at face value, our Y p 
implies a baryon-to-photon number ratio r) = (4.7t! 0 ° 8 )x 10” 10 and a baryon mass fraction IWiioo 
= 0.017±0.005 (2 <t), consistent with the values obtained from deuterium and Cosmic Microwave 
Background measurements. Correcting f p upward by 2^t% would make the agreement even better. 
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1. Standard big bang nucleosynthesis 

The standard hot big bang model of nucleosynthesis (SBBN) is one of the key 
quantitative tests of big bang cosmology, along with the Hubble expansion and 
the cosmic microwave background radiation. In the SBBN, four light isotopes, D, 
3 He, 4 He and 7 Li, were produced by nuclear reactions a few minutes after the birth 
of the Universe. Given the number of relativistic neutrino species and the neu¬ 
tron lifetime, the abundances of these light elements depend on one cosmological 
parameter only, the baryon-to-photon ratio r/, which in turn is directly related to 
the density of ordinary baryonic matter fij,. The ratio of any two primordial abun¬ 
dances, for example that of D to H gives r/, and accurate measurements of the other 
three light elements, for example 4 He/H, tests SBBN. 

Of all light elements, the abundance of deuterium (D) is the most sensitive to 
the baryonic density. The primordial D abundance can be measured directly in 
low-metallicity absorption line systems in the spectra of high-redshift quasars. The 
quasar is used as a background light source, and the nearly primordial gas doing 
the absorbing is in the outer regions of intervening galaxies or in the intergalactic 
medium (the so-called Lyman a clouds). Tytler and his group (see Tytler et al. 
2000 for a review) have vigorously pursued this type of measurements. They have 


2 


THUAN & IZOTOV 


now obtained D/H measurements in the line of sight towards 4 quasars. Combining 
all measurements, they found all their data arc consistent with a single primordial 
value of the D/H ratio: (D/H) p = 3.0±0.4x 10~ 5 (O'Meara et al. 2001). This latest 
value is about 10% lower than their previous value (D/H) p = 3.39±0.25xl0 5 
(Buries & Tytler 1998). 

The primordial abundance of 3 He is also quite sensitive to the baryon density, 
though less than the D abundance. It has not been yet measured, mainly because 
low-mass stars make a lot of 3 He, increasing its value in the interstellar medium 
of the Milky Way well above the primordial value. Furthermore, the amount of 
3 He destroyed in stars is unknown. Bania, Rood & Balser (2000) have measured 
an average 3 He/H = 1.3±0.6x10 3 in Galactic H II regions. This value represents 
the average in the interstellar medium of the Milky Way, but there exists no good 
way to extrapolate the 3 He abundance to the primordial value. 

Old halo stars that formed from nearly pristine gas with very low iron abun¬ 
dances during the gravitational collapse of the Milky Way show approximate¬ 
ly constant 7 Li/H (the so called “Spite plateau”, Spite & Spite 1982), implying 
that their 7 Li is nearly primordial. Creation or depletion of 7 Li may make the 7 Li 
abundances of halo stars deviate from the primordial value. Creation of 7 Li in the 
interstellar medium by cosmic ray spallation prior to the formation of the Milky 
Way has to be less than 10-20%, so as not to produce more Be than is observed 
(Ryan, Norris & Beers 1999). There is still considerable debate concerning the 
possible depletion of 7 Li at the surface of stars. Depletion mechanisms that have 
been proposed include mixing due to rotation or gravity waves, mass loss in stellar 
winds and gravitational settling. Depending on the exact depletion mechanism, the 
primordial lithium abundance varies from ( 7 Li/H) p = (1.73±0.21)x 10~ 1(l (Boni¬ 
facio & Mol am 1997) to (2.24±0.57)x 10 10 (Vauclair & Charbonnel 1998), to 
(3.9±0.85)x 10 -10 (Pinsonneault et al. 1999). 

Because of the relatively large uncertainties in the determination of the primor¬ 
dial abundances of 3 He and 7 Li, the primordial abundance of 4 He plays a key role 
for deriving O/, independently of D measurements, and is crucial for checking the 
consistency of SBBN. We discuss next how the primordial 4 He mass fraction T p is 
determined, and the uncertainties which enter in such a determination. 


2. The primordial 4 He abundance as derived from blue compact dwarf 

galaxies 

Because of the relative insensitivity of 4 He production to the baryon ic density of 
matter, T p needs to be determined to a precision of about one percent to provide 
useful cosmological constraints. This precision can in principle be achieved by 
obtaining high signal-to-noise ratio spectra of a class of star-forming dwarf galax¬ 
ies called Blue Compact Dwarf (BCD) galaxies. These are low-luminosity (Mb 
> -18) systems undergoing an intense burst of star formation in a very compact 
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region (less than 1 kpc) which dominates the light of the galaxy and which shows 
blue colors and a H II region-like emission-line optical spectrum (Thuan & Izo¬ 
tov 1998a). BCDs arc ideal laboratories in which to measure the primordial 4 He 
abundance because, with an oxygen abundance O/H ranging between 1/50 and 1/3 
that of the Sun, BCDs arc among the most metal-deficient gas-rich galaxies known. 
Their gas has not been processed through many generations of stars, and thus best 
approximates the pristine primordial gas. Izotov & Thuan (1999) have argued that 
BCDs with O/H less than ~ 1/20 that of the Sun may be genuine young galax¬ 
ies, with ages less than several 100 Myr. Thus Y p can be derived accurately in very 
metal-deficient BCDs with only a small correction for Helium made in stars. More¬ 
over, the theory of nebular emission is well understood enough to allow to convert 
He emission-line strengths into abundances with the desired accuracy. 



4000 ' ' ' ' 5000 ' ' ' 6000 ' ' 7000 

Wavelength (A) 


Figure 1. MMT spectra of the NW (top) and SE (bottom) regions of I Zw 18. It is evident that 
underlying He I stellar absorption is much more important in the NW than in the SE component. All 
marked He I lines in the spectrum of the SE component are in emission, while the He I 4026 and 
4921 lines are in absorption and the He I 4471 emission line is barely seen in the spectrum of the 
NW component. 


Tp is generally determined by linear extrapolation of the correlations Y- O/H 
and Y-N/H to O/H = N/H = 0 as first proposed by Peimbert & Torres-Peimbert 
(1974), where Y, N/H and O/H arc respectively the 4 He mass fraction, the Oxygen 
and Nitrogen abundances relative to Hydrogen of a sample of dwarf irregular and 
BCD galaxies. Based on a relatively large sample of ~ 45 BCDs, we (Izotov, Thuan 
& Lipovetsky 1994, 1997, hereafter ITL94, ITL 97; Izotov & Thuan 1998, here¬ 
after IT98) have obtained Y p = 0.2443+0.0015 with dY/dZ = 2.4+1.0 (see Thuan 
& Izotov 1998b, 2000 for reviews). This result is quite robust as it is in excellent 
agreement with the average Y p = 0.2452+0.0015 (Izotov et al. 1999) determined in 
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the two most metal-deficient BCDs known, I Zw 18 (Zq/50) and SBS 0335-052 
(Zq/ 41), where the correction for He production is smallest. The quoted statistical 
error (la) is ^ 1% because of the high quality of the spectra and the large size of 
the sample. 

Our derived value of Y p is significantly higher than those of 0.228±0.005 by 
Pagel et al. (1992), of 0.234±0.002 obtained by Olive, Steigman & Skillman (1997), 
of 0.2345T0.0026 by Peimbert, Peimbert & Ruiz (2000) and of 0.2384T0.0025 
by Peimbert, Peimbert & Luridiana (2001). These values are lower than ours by 
several times the quoted statistical error and lead to very different cosmological 
consequences. This implies that systematic effects in the determination of V p play 
an important role and need to be considered carefully to achieve the desired accu¬ 
racy. We discuss the most important systematic effects below, specifically in the 
two most metal-deficient BCDs known, I Zw 18 and SBS 0335-052. Because dif¬ 
ferent systematic effects arc at play in these two objects, they allow to illustrate the 
relative importance of each nicely. 


3. Systematic effects 


3.1. Reddening 

One possible systematic effect may be due to the adopted interstellar extinction 
curve. We use the Galactic extinction curve by Whitford (1958). There is evidence 
that the extinction curve does change when metallicity decreases, being steeper at 
short wavelengths (e.g. Rocca-Volmerange et al. 1981), but the changes are mainly 
in the ultraviolet and arc small in the optical range. Relying mainly on the three 
Baimer line ratios Hcr/H/3, Hy/H't and Hh/H/t, we have adopted an iterative pro¬ 
cedure to derive simultaneously both the extinction and the absorption equivalent 
width for the hydrogen lines assumed to be the same for all lines (ITL94). Since 
the derived extinction correction is applied to the He I lines, any uncertainty in it 
will propagate in the derivation of the final He abundance. A stringent observation¬ 
al check for the adequacy of the extinction curve is the good agreement between 
the corrected intensities of the Balmer hydrogen emission lines and the theoretical 
values for hydrogen recombination line intensities. Using Monte-Carlo simulations 
of the hydrogen Balmer ratios, Olive & Skillman (2001) estimate the uncertainties 
in the reddening corrections to be about 1-2 % for the blue He lines and 3-4 % for 
the red ones. However, the uncertainties for the red lines can be substantially less 
if the corrected Ho/H/7 ratio matches the theoretical ratio. 

3.2. Underlying stellar absorption in He i lines 

Underlying stellar absorption in He I lines caused by hot stars can decrease the 
intensities of nebular He I lines. Model calculations of synthetic absorption line 
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Figure 2. Spatial distributions of the He I nebular emission line equivalent widths in I Zw 18 (left 
panel) and in SBS 0335-052 (right panel). The error bars are la deviations. The value of the mini¬ 
mum equivalent width for each He I emission line is given. 


strengths in star forming regions by Olofsson (1995) show that the equivalent 
widths of He I absorption lines decrease as the starburst ages. Furthermore, the 
dependence of He I equivalent widths on metallicity is small, and the equivalent 
width of the He I 4471 absorption line can be as high as 0.35A. Unfortunately, 
similar calculations for other important lines used in the determination of Ip such 
as He I 5876, 6678 and 7065 are not yet available. The effect of underlying He I 
stellar absorption is most important for the emission lines with the smallest equiva¬ 
lent widths. Therefore, the He I 5876 emission line which has the largest equivalent 
width is the least affected by such absorption, while the He I 4471 emission line 
is the most affected because of its low equivalent width, the effect of underlying 
absorption being 5-10 times larger for it than for the 5876 line. 

Figure 1 shows the spectra of the 2 brightest centers of star formation in I Zw 
18, the NW and SE components, with the He I lines marked. It is clear that the 
NW component suffers far more stellar absoiption than the SE one: all marked He 
I lines of the SE component arc in emission while the two He 14026 and 4921 lines 
arc in absorption and the He I 4471 line is hardly detected in the NW component, 
its intensity being decreased by a factor of ~ 2 (Izotov et al. 1999). However, the 
strong He I absorption in the NW component of I Zw 18 constitutes the exception 
rather than the rule. It is relatively less important in other very metal-deficient 
BCDs, such as SBS 0335-052, because the equivalent widths of their He I emission 
lines arc considerably larger than those of the NW component of I Zw 18. This can 
be seen in Figure 2 which compares the spatial distributions of the He I nebular 
emission line equivalent widths (EW) in I Zw 18 and SBS 0335-052. While the 
maximum EW s of the He I emission lines in the SE component of I Zw 18 arc 


thuan.tex; 1/02/2008; 20:59; no v.; p.5 




6 


THUAN & IZOTOV 


close to those in the central brightest part of SBS 0335-052, the He I line EW s 
in the NW component of I Zw 18 arc several times smaller. The largest ratio of 
minimum values of equivalent widths in I Zw 18 as compared to in SBS 0335-052 
is ~ 14 for the 4471 line. The effect of underlying stellar absorption is smaller 
for the other He I lines because of their higher equivalent widths. Thus for the 
6678 line, Izotov et al. (1999) estimate it to be ~ 5% in the NW component of 
I Zw 18, but to be less than 1% in its SE component and in SBS 0335-052. As 
for the 5876 line, it is less than 0.4% in SBS 0335-052 (it is contaminated by the 
Galactic Na I 5890 line in I Zw 18, so is unusable). It was mainly because the 
importance of underlying He I stellar absorption was insufficiently recognized in 
the NW component of I Zw 18 that led to the low Tp values of Pagel et al. (1992) 
and Olive et al. (1997). 



SW<->NE • He I X5876 

□ He I A6678 


b) N (He + ) = N (S + ) SBS 0335-052 



Distance along the slit (arcsec) 


Figure 3. The spatial distributions of the helium mass fractions in SBS 0335-052 derived from the 
He I 4471, 5876 and 6678 emission line intensities. The intensities of the He I emission lines in a) 
are corrected for fluorescent and collisional enhancement with an electron number density 7V e (He 
II) and an optical depth r(3889) derived self-consistently from the observed He I 3889. 4471. 5876. 
6678 and 7065 emission line intensities. The points for the 4471 line are below the points for the 
other lines because of underlying He I stellar absorption. For comparison, the intensities of the He I 
emission lines in b) are corrected only for collisional enhancement, with an electron number density 
N e (S II). The points for the different lines do not agree anymore. The la error bars are shown only 
for the He mass fraction derived from the He I A5876 emission line. They are larger in b) because of 
the large uncertainties in the determination of N e (S II). 


In our work thus far, we have not corrected self-consistently for underlying He 
I self-absorption. We have simply not used regions where underlying absorption is 
important (such as the NW component in I Zw 18), or not averaged in lines that give 
a Y clearly deviant from the Y s from other lines. Thus, in the case of SBS 0335- 
052, Figure 3a shows that the Y obtained from He I 4471 line is systematically 
below the values from the He I 5876 and 6678 lines. It is affected by underlying 
absorption and is not included. In the future, we plan to solve for He I absorption 
(assumed to be the same for all He I lines) self-consistently, as we did for the H I 
lines (see below). 
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3.3. HE I AND H I EMISSIVITIES 

Line ratios corrected for reddening and absoiption are converted to He/H abun¬ 
dance ratios by using theoretical emissivities calculated from recombination and 
radiative cascade theory. We use the H I emissivities of Brocklehurst (1971) which 
arc in excellent agreement with the more recent ones by Hummer & Storey (1987) 
for the range of temperatures and densities in metal-deficient BCDs, and the He I 
emissivities of Smits (1996). Benjamin, Skillman & Smits (1999) have estimated 
that uncertainties in the theoretical He I emissivities can be as large as 1.5%, ~ 3 
times worse than the accuracy expected from comparing Brocklehurst (1972) and 
Smits (1996) emissivities. 

3.4. Electron density determination 

To determine element abundances, we adopt a two-zone photoionized H II region 
model (ITL94,ITL97,IT98): a high-ionization zone with temperature T e (O III), 
where the O III, Ne III and Ar IV lines originate, and a low-ionization zone with 
temperature T e (0 II), where the O II, N II, S II and Fe III lines originate. As for the 
Ar III and S III lines they originate in the intermediate zone between the high and 
low-ionization regions. The [S II] A6717/A6731 line ratio is used to determine the 
electron density N e (S II) according to five-level atom calculations. 

Previous authors have set A). (He n) = N e (S II). However, this is not appropriate 
as N e (S II) measures the density in the low-ionization zone, while He II is produced 
in the high-ionization zone. Furthermore, the [S II] ratio is fairly insensitive for 
densities below 100 cm~ 3 . Thus it is much better to determine A). (He n) directly 
from the He I lines themselves, in a self-consistent manner, which we have done. 
Because the electron density enters linearily in the calculation of the correction 
for collisional enhancement of the He I emission lines (see below), its different 
estimate by us and previous authors is also responsible for the difference between 
our Tp and theirs. 

3.5. The collisional excitation of hydrogen lines 

It is also generally assumed that case B recombination theory holds, i.e that the line 
photons (usually resonance lines of abundant ions) are scattered so many times that 
their downward radiative transitions can effectively be neglected. In other words, 
the H II regions are considered to be optically thick in the Lyman transitions but 
optically thin in the other transitions of both hydrogen and helium atoms. (The 
opposite assumption, case A, where all emitted photons escape without absorption 
and there is no radiative transfer problem, does not apply to nebulae with large 
enough amounts of gas and optical depths to be seen). However, there arc known 
physical effects that make the helium lines deviate from case B. One such effect is 
the collisional excitation of the B aimer hydrogen lines by thermal electrons. 
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Figure 4. The spatial distributions of the He I and He II 4686 nebular emission line intensities in I 
Zw 18 (left panel) and in SBS 0335-052 (right panel). The error bars are la deviations. 


This process and its effect on the determination of Ip was first discussed by 
Davidson and Kinman (1985). Low-metallicity H II regions have high enough elec¬ 
tron temperatures so that the observed fluxes of lines like Ha may be overestimated 
by this effect, leading to an underestimate of Y p . Davidson & Kinman (1985) esti¬ 
mated crudely this effect to be ~ 2% for I Zw 18. Sasselov & Goldwirth (1995) 
examined this effect with detailed radiative transfer calculations and also found it 
to cause an increase of ~ 2% in the Y of I Zw 18 and up to 3% for other metal- 
deficient BCDs. Stasiriska & Izotov (2001) used a grid of photoionization models 
to show that the effect of collisional excitation on the Ha/H/3 ratio can be as high as 
8%, resulting in an upward correction in Y of up to 5% in objects like SBS 0335- 
052. Thus this effect can be one of the most important sources of systematics in the 
determination of Tp. 

3.6. The collisional and lluorescent enhancements of He i lines 

In the high range of electron temperatures found in metal-deficient BCDs (10000- 
20000 K), collisional excitation from the metastable level 2 3 S level of He I can 
be important in populating the higher levels and making the He I line intensities 
deviate from pure recombination values. Another effect that also leads to deviation 
from pure recombination values is self-absorption in some optically thick emission 
lines which populates the upper levels of He I, a mechanism called fluorescence. 
The emission lines most sensitive to fluorescence in the optical range are the He I 
3889 and 7065 lines. The He I 7065 also plays an important role because it is par¬ 
ticularly density sensitive. In contrast to collisional excitation, which increases the 
line intensities of all He I lines, the fluorescent mechanism decreases the intensity 
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of the He I 3889 line as its optical depth increases, while increasing the intensities 
of the He I 4471, 5876, 6678, and 7065 lines. 



o L—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i— 
5500 6000 6500 7000 
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Figure 5. Keck spectrum of the BCD SBS 0335-052 with labeled He I emission lines. 

How do these two effects play out in I Zw 18 and SBS 0335-052? They are 
less important in I Zw 18 than in SBS 0335-052 because the physical conditions 
in the two BCDs arc quite different. While in I Zw 18 the electron number density 
is small {N e (S II) < 100 cm ’) and collisional enhancement has a minor effect on 
the derived helium abundance, the electron number density in SBS 0335-052 is 
considerably higher ( N e (S II) ~ 500 cnr 3 in the central part of the H II region). 
Additionally, the linear size of the H II region in SBS 0335-052 is ~ 5 times larger 
than in I Zw 18, making it optically thick for some He I transitions. Therefore, both 
collisional and fluorescent enhancements of He I emission lines play a significant 
role in this galaxy (Izotov et al. 1999). Figure 4 shows that the spatial distribution 
of He I emission line intensities in SBS 0335-052 is very different from that in I 
Zw 18. The increase of He I 5876 and He I 7065 emission line strengths by ~ 20% 
and ~ 75% respectively in the central part of SBS 0335-052 within a radius ~ 
2" is caused by collisional and fluorescent enhancement. The increase of the He I 
A6678 emission line intensity is only < 4%. The combined effect of collisional 
enhancement and underlying stellar absorption results in a small depression in the 
He I 4471 intensity in the central region. As for I Zw 18, the main effect is under¬ 
lying He I stellar absorption causing the dip in the He I emission line intensities at 
the location of the NW component. 

To correct the He I line intensities for collisional enhancement, we use the cor¬ 
rection factors calculated by Kingdon & Ferland (1995) based on collisional rates 
by Sawey & Berrington (1993). To correct the He I line intensities for fluorescent 
enhancement, we have fitted the Robbins (1968) correction factors with polyno¬ 
mials as given in IT98. Since the collisional enhancement factor of the He I lines 
depend exponentially on the electron temperature and linearly on the electron den¬ 
sity N e (He II), we correct for these two effects and determine N e (He II) at the 
same time, in a self-consistent way so that the He I 5876/4471, 6678/4471, and 
7065/4471 line ratios have their recombination values, after correction for both 
collisional and fluorescent enhancement. 
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3.7. The non-coincidence of the H + and He + Stromgren spheres 

Depending on the hardness of the ionizing radiation, the radius of the He + sphere 
can be smaller or larger than the radius of the H + sphere. When the ionizing radia¬ 
tion is soft, the first case prevails, it is necessary to consider the presence of unseen 
neutral helium in the H + zone and a correction needs to be made, resulting in a 
higher He abundance. On the other hand, if the ionizing radiation is hard, the sec¬ 
ond case holds, there is neutral hydrogen in the He + zone, which results in a down¬ 
ward correction of the He abundance. The correction in Y can be as high as several 
percent, either upward or downward (e.g. Steigman, Viegas & Gruenwald 1997, 
Peimbert, Peimbert & Ruiz 2000, Sauer & Jedamzik 2001). The hardness of the 
radiation is usually characterized by the “radiation softness parameter” rj defined 
by Vilchez & Pagel (1988) as r/ = §r^++, and the ionization parameter U. For both 
I Zw 18 and SBS 0335-052, the ionizing radiation is hard, and the correction to Y 
is downward. Using the extensive grid of correction factors as functions of 7] and 
U, calculated by Sauer & Jedamzik (2001) using photoionized H II region models, 
we found the downward correction to be less than 1 % for both BCDs. 

3.8. Temperature structure 

To convert the He I line intensities into abundances, we have set the electron tem¬ 
perature T e (He II) to be equal to T e ( O III) as obtained from the [O III] 4363 / 
(4959+5007) ratio. However as emphasized by Peimbert, Peimbert & Ruiz (2000), 
detailed modeling of BCDs like I Zw 18 (Stasinska & Schaerer 1999) and examina¬ 
tion of photoionization models (Stasinska 1990) suggest that T e (He II) is smaller 
than T e (O III) by at least 5% in this type of object (deviations from a constant 
temperature are sometime called “temperature fluctuations”). Peimbert, Peimbert 
& Luridiana (2001) found that this leads to a downward correction of Y of about 
3%. However, examination of the correction factors for temperature effects derived 
by Sauer & Jedamzik (2001) using photoionization models give dow nward correc¬ 
tions of Y that arc considerably smaller for I Zw 18 and SBS 0335-052, less than 
1 %. 


4. Cosmological implications 

In summary, the important systematic effects that may affect our determination of 
Yp from spectra of the two most metal-deficient BCDs known, I Zw 18 (the SE 
component) and SBS 0335-052, and that have not been taken into account by our 
self-consistent procedure arc: 1) the collisional excitation of Hydrogen lines that 
can increase Y p by up to 5%; 2) the non-coincidence of the H + and He + zones that 
may decrease Y p by ~ 1%; 3) the temperature fluctuations that may decrease it by 
~ 1 - 3%; and 4) underlying stellar He I absorption that may increase it by ~ 1%. 
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Thus, combining all those systematic effects, our Ip value may be underestimated 
by as much as 2 - 4%. 

Taken at face value, our >' p = 0.2452±0.0015 implies a baryon-to-photon num¬ 
ber ratio // = (4. 7 -o°8)x 10 10 . This translates to a baryon mass fraction Qi ) h 2 ]()() = 
0.017±0.005 (2a) where h\oo is the Hubble constant in units of 100 km s _1 Mpc ' 1 . 
This value is consistent with the one of 0.020T0.002 (2a) derived from the primor¬ 
dial deuterium abundance measured in high-redshift hydrogen clouds backlit by 
distant quasars (O’Meara et al. 2001), vindicating SBBN. It is also consistent with 
the baryon mass fraction Qi,h 2 ]{)() = 0.022±0.003 (I a) inferred from measurements 
of the angular power spectrum of the Cosmic Microwave Background (CMB, Net- 
terfield et al. 2001). Note that correcting upward Y' p by 2M% would bring it into 
even better agreement with the deuterium and CMB measurements. 


5. Future work 

It is cleai - from the previous discussion that, to decrease or eliminate the main 
systematic effects in the determination of Y p , it is best to determine all the four 
following quantities - the electron density N e (He II) and temperature T e (He II) of 
the He II zone, the optical depth r (3889) in the He I 3889 line, and the equivalent 
width for underlying He I absorption - in a totally self-consistent manner. ITL94, 
ITL97 and IT98 have used the five He I 3889, 4471, 5876, 6678 and 7065 lines 
to solve self-consistently for N e (He II) and r (3889). We plan to add the He I 
4026, 4438 and 4922 lines, giving a total of 8 lines, to solve self-consistently for 
all four quantities. The last three lines are particularly sensitive to underlying He I 
stellar absorption. However their intensities arc less than 1/10 that of the 5876 line, 
and very high signal-to-noise ratio spectra arc needed to determine their intensities 
precisely. Such a spectrum obtained with the Keck telescope and where all 8 He I 
lines arc marked is shown in Figure 4 for the BCD SBS 0335-052. 

Future work to improve on the determination of Y p will consist of: 1) obtain¬ 
ing deep spectra for the most metal-deficient BCDs known, to determine Y p using 
a self-consistent method based on 8 lines; 2) detailed modeling of BCDs such as 
done for I Zw 18 by Stasinska & Schaerer (1999). It was found that a simple pho¬ 
toionization model is insufficient to account for the high [O III] 4363/5007 ratio 
in I Zw 18 and that an additional heating mechanism such as shock heating is 
necessary. Such a modeling will also help to quantify potential temperature fluctu¬ 
ations; and 3) searching for more extremely metal-deficient BCDs such as I Zw 18 
and SBS 0335-052, to increase the number of objects where we can determine Y p 
without a large correction for He production by stars. 
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